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I  INTRODUCTION 


The  purpose  of  this  report  is  to  describe  a  technique  for  solving  steady  state 
coupled  structure-acoustics  problems  that  has  been  implemented  by  the  transducer 
modeling  group  (Code  712)  at  the  Naval  Ocean  Systems  Center  (NOSC).  This  tech¬ 
nique  is  applicable  to  both  structural  radiation  problems  and  structural  scattering 
problems. 

The  structural  portion  of  the  problem  is  modeled  using  a  standard  finite  element 
program.  The  particular  program  used  at  NOSC  is  MARTSAM,  but  virtually  any 
general  finite  element  code  could  be  used.  The  acoustical  portion  of  the  problem  is 
modeled  using  the  Combined  Helmholtz  Integral  Equation  Formulation  (CHIEF)1. 
This  technique  has  been  implemented  at  NOSC  in  a  computer  program  that  is  also 
called  CHIEF2.  The  original  documentation  on  CHIEF  described  only  its  acoustic 
radiation  capabilities,  but  the  capability  to  do  scattering  problems  has  been  imple¬ 
mented  and  is  described  in  the  revised  CHIEF  Users  Manual.  The  coupling  of  the 
structural  and  acoustic  models  is  performed  in  a  separate  program  that  combines 
certain  outputs  of  the  finite  element  and  CHIEF  programs.  The  approach  that  has 
been  implemented  follows  that  described  in  Wilton3  except  for  the  way  in  which  the 
continuity  of  normal  velocity  is  enforced. 

This  report  is  organized  as  follows:  Section  II  summarizes  the  basic  finite  element 
structural  relations;  Section  III  summarizes  the  basic  CHIEF  acoustic  relations;  Sec¬ 
tion  IV  describes  the  combination  of  the  CHIEF  and  finite  element  results  to  solve 
the  coupled  problem;  Section  V  describes  the  output  quantities  that  are  computed  for 
a  scattering  problem,  and  compares  numerical  results  for  scattering  from  an  elastic 
spherical  shell  to  the  analytical  solution;  and  Section  VI  describes  the  output  quan¬ 
tities  that  are  computed  in  a  radiation  problem,  and  compares  numerical  results  for 
radiation  from  a  locally  driven  elastic  spherical  shell  to  the  analytic  series  solution. 
Section  VII  contains  a  summary  of  the  features  of  the  model  described  in  this  report. 

*H.  Schenck,  Improved  Integral  Formulation  for  Acoustic  Radiation  Problems,  J.  Acoust.  Soc. 
Am.  44,  41-58  (1968). 

2G.W.  Benthien,  D.  Barach,  and  D.  Gillette,  CHIEF  Users  Manual,  Naval  Ocean  Systems 
Center  Technical  Document  970,  Rev.  1,  November  1988. 

3D.T.  Wilton,  Acoustic  Radiation  and  Scattering  from  Elastic  Structures,  International  Journal 
for  Numerical  Methods  in  Engineering  13,  123-138  (1978). 
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II  STRUCTURAL  EQUATIONS 


The  elastic  structure  occupies  a  region  R  that  is  bounded  by  a  surface  dR. 
The  vibration  of  the  structure  is  assumed  to  be  described  by  the  equations  of  linear 
elasticity.  Since  the  equations  are  linear,  the  time  dependence  of  all  variables  will  be 
expressed  in  terms  of  Fourier  transforms,  i.e., 

u(x,0  =  /  U(z,u>)e,u"du;, 

where  the  complex  function  u(z,u>)  is  called  the  Fourier  transform  of  u(z,  t).  The 
linearity  of  the  equations  of  motion  allows  them  to  be  expressed  entirely  in  terms  of 
the  complex  Fourier  components.  In  what  follows,  the  explicit  dependence  of  variables 
on  the  frequency  uj  will  usually  be  omitted. 

In  the  finite  element  method,  summarized  here,  the  displacement  u  is  approxi¬ 
mated  by  a  linear  combination  of  vector  interpolation  functions  d>m(z),  i.e., 

u(z)  =  Y,  Um<t>m{x)  ,  *  €  R.  (1) 

m=l 

Typically,  the  values  Um  correspond  to  the  components  of  displacement  at  a  finite 
number  of  points  in  the  body  called  nodes,  and  the  functions  <f>m{x)  form  a  basis  for 
piecewise  polynomial  interpolation  between  nodes.  The  finite  element  equations  can 
be  obtained  by  utilizing  the  approximation  of  equation  (1)  in  a  variational  formulation 
of  the  equations  of  elasticity4.  The  resulting  finite  element  equations  are  usually 
expressed  in  the  matrix  form 

(~u2M  +  K)U  =  Fl,  (2) 

where  the  mass  matrix  M  has  the  components 

M*  =  /  P*<t>,  ■  <t>,  dV  (3) 

J  H 

p, . . .  mass  density  of  material; 

4O  C.  Zienkiewicz,  The  Finite  Element  Method,  Third  Edition,  McGraw-Hill,  1977,  Chap¬ 
ter  3. 
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the  stiffness  matrix  K  has  the  components 


Kt}  =  /  C,[V<t>,)-<7<t>:dV  (4) 

J  H 

V<t>  =  \{V<t>  +  V<t>T) 

C,  ...elasticity  tensor  relating  the  stress  tensor  and  strain  tensor; 
and  the  load  vector  FL  has  the  components 


t . . .  stress  vector  load  on  dR. 

The  unknown  vector  V  consists  of  all  displacement  components  at  all  the  nodes  in  R. 


For  a  structure  immersed  in  a  fluid  medium  it  is  convenient  to  write  the  load 
vector  FL  as 

Fl  =  FR  +  Fd ,  (6) 

where  FR  is  associated  with  the  acoustic  loading  of  the  fluid  medium  and  F°  is 
associated  with  the  other  known  forces  which  drive  the  structure. 

Equations  (2)  and  (6)  can  be  combined  to  yield 

(~v2M  +  K)U  =  FR  ■+  Fd.  (7) 

For  an  ideal,  nonviscous  fluid  medium,  the  stress  vector  t  can  be  written  in  terms  of 
the  acoustic  pressure,  i.e.. 

t  =  -pn  (8) 

n . . .  unit  outward  normal  vector  (pointing  into  fluid). 


If  S  is  the  portion  of  dR  that  is  in  contact  with  the  fluid,  then  the  load  vector  FR 
can  be  written  as  follows: 

F*  =  -  JsP4>m-ndS.  (9) 

In  Section  III,  the  governing  equations  for  the  fluid  medium  will  be  discussed.  These 
equations  will  relate  the  pressure  p  in  equation  (9)  to  the  normal  velocity  on  the 
surface  5. 
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Ill  ACOUSTIC  EQUATIONS 


In  a  scattering  problem  one  is  interested  in  the  scattered  pressure  field  p,  due 
to  a  plane  wave  p,nc  impinging  on  the  surface  S  bounding  the  elastic  structure.  The 
total  pressure  p  is  related  to  the  scattered  pressure  p,  by 

P  =  Pinc  +  P.-  (10) 


In  a  radiation  problem  there  is  no  incident  wave  p,nc,  so  p  and  p,  are  the  same.  It  is 
shown  by  Wilton5  that  the  total  pressure  p  satisfies  the  Helmholtz  integral  formulas 


dGjry) 

dnr 


dpjr) 

dnr 


G(rS) 


dS(r) 


jTP(r')  -  Pine(r'), 

4  ir 


r'  on  5 


(  Ha) 


—  Pmc(r  )> 

r'  interior  to  S  (lib) 


=  P(r')  -  pinc(r'), 

r'  external  to  5 


(  He) 


where 


G(r,r')  = 


pi*|r-r'| 


4ir|r  —  r'l 

is  the  three-dimensionai  free  Green’s  function  for  the  Helmholtz  differential  equation 


V2p  +  Jt2p  -  0; 

a  is  the  solid  angle  at  r1  £  S  subtended  by  the  fluid  exterior  to  5;  and  ^  denotes 
differentiation  in  the  outward  normal  direction  n  at  r. 


At  any  point  r  6  S  at  which  the  normal  n(r)  is  well-defined,  the  normal  derivative 
of  the  pressure  p  is  related  to  the  normal  velocity  v(r)  by6 

=  -iwpn(r),  (12) 

5D.T.  Wilton,  ibid,  page  125. 

6D.T.  Wilton,  ibid,  page  126. 
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where  p  is  the  fluid  density.  Combination  of  equations  (11)  and  (12)  gives 

^P(r')  -  Jsp(r)d~£r)  dS{  r)  =  iup  Jsv(r)G{r,r')  dS(r)  +  pinc{r')  (13a) 
for  r'  on  the  surface  5; 

-  dS(r)  =  i^P  Jsv{r)G{r,r')  dS{r)  +  pinc{r')  (13b) 

for  r'  interior  to  the  surface  5;  and 

P,(r  )  =  Is  {P(r)~in  ^  +  iwMr)G(r,r')  j  d5(r)  (13c) 

for  r'  exterior  to  the  surface  5. 

If  v  and  pinc  are  prescribed  on  5,  then  equation  (13a)  becomes  an  integral  equation 
of  the  second  kind  for  p  on  S.  This  integral  equation  is  well  suited  for  numerical 
computations  for  most  wavenumbers  k ,  but  its  solution  is  not  unique  whenever  k  is  an 
eigenvalue  of  the  interior  Dirichlet  problem.  One  effective  approach  to  overcoming  this 
lack  of  uniqueness  is  to  use  the  combined  Helmholtz  integral  equation  formulation7. 
In  this  approach,  an  overdeter  mined  system  of  equations  is  formed  by  supplementing 
the  integral  equation  (13a)  with  the  integral  equation  (13b)  evaluated  at  selected 
points  r'  interior  to  S. 

The  combination  of  equations  (13a)  and  (13b)  can  be  solved  numerically  by  ap¬ 
proximating  p  and  u  on  5  by  piecewise  constant  functions,  i.e., 

p(r)  =  pn  and  v{r)  =  vn  for  r  €  5„,  (14) 

where  Si,...,Sjv  are  nonoverlapping  subregions  of  S  whose  union  is  S.  Collocation 
of  (13a)  at  the  N  surface  points  rm  (m  =  1  ,...,/V)  where  the  surface  is  smooth 
(a  =  2tt),  and  (13b)  at  the  N'  interior  points  rm  (m  =  N  +  1, . . . ,  N  +  N')  leads  to 
the  matrix  relation 

AP  =  BV  +  Ptnc,  (15) 

'  H.  Schenck,  ibid. 
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where 


Amn  =  6mn-f  dG(r’^l  dS(r)  m  =  \,...,N  +  N'  n  =  l . TV  (16) 

is„  onr 

J5mn  =  iu>p  f  G(r,  rm)  dS(r )  m  =  1, . . . ,  N  +  TV'  n  =  l,...,N  (17) 

•/5n 


!  Pi  ^ 

^  Ui  ^ 

1  Pxnc(r  l)  N 

p  = 

j 

v  = 

| 

P  — 

»  1  tnc 

| 

\  PN  ) 

l  VN  j 

\  pinc{rN+N')  ) 

Equations  15-17  represent  a  special  case  of  more  general  relations  given  by  Wilton8. 
In  the  next  section,  the  acoustic  relation  (15)  will  be  coupled  with  the  finite  element 
structural  equations.  Solution  of  the  coupled  equations  will  determine  the  surface 
pressures  jq, . . .  ,pN  and  the  normal  velocities  u,, . . . ,  v /y. 


Applying  the  piecewise  constant  approximation  to  equation  (13c)  results  in 

P«(r')  =  ]CP"/  dG(r'.r  )  dS(r)  +  ivp  ^2  vn  f  G(r,  r1)  dS(r).  (18) 

n  =  l  JS n  OTIt  n=l  J  S* 


Equation  (18)  can  be  used  to  calculate  the  scattered  pressure  p„  at  any  point  r' 
exterior  to  S.  If  is  an  arbitrary  set  of  field  points  where  it  is  desired  to 

calculate  the  scattered  pressure,  then  equation  (18)  reduces  to  the  matrix  equation 


P.  =  afp  +  bfv , 

(19) 

where 

(20) 

and 

Bmn  =  '“P  f  G{r,rm)dS(r). 

Js„ 

(21) 

SD.T.  Wilton,  ibid,  pages  126-127. 
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Figure  1.  Far-field  quantities. 


Often  it  is  desired  to  calculate  the  scattered  or  the  radiated  pressure  (Pinc  —  0)  in  the 
‘far-field’  of  S  where  the  field  pressure  spreads  spherically.  Let  us  define  a  normalized 
far-field  pressure  pFF  by 


(22) 


where  R  is  the  distance  from  an  origin  near  the  surface  S  to  the  field  point.  It  can 
be  shown9  that  equation  (19)  can  be  approximated  in  the  far-field  by 


pFF  _  AFFp  +  BFFy • 


(23) 


where 


and 


i  FF 


-£/<*. 

47T  Jsn 


•  n)e 


iiRmi(r) 


dS(r) 


B 


FF 

mn 


4ir  Jsn 


(24) 

(25) 


The  vectors  Rm  and  6(r)  are  pictured  in  figure  1  and  Rm  is  the  unit  vector  in  the 
direction  Rm.  An  important  point  to  note  is  that  the  matrices  AFF  and  BFF  are 
functions  of  6  and  <f>,  but  are  independent  of  R. 


9G.W.  Benthien,  D.  Barach,  and  D.  Gillette,  ibid,  Appendix  A4. 
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The  acoustic  relations  presented  in  this  section  have  been  incorporated  into  a 
computer  program  called  CHIEF.  This  program  was  initially  set  up  to  handle  only 
acoustic  radiation  problems,  but  it  has  recently  been  modified  to  handle  scattering 
problems  as  well.  The  CHIEF  program  has  special  provisions  to  take  advantage  of 
one,  two,  or  three  planes  of  symmetry  and  finite  degrees  of  rotational  symmetry  in 
the  geometry  and/or  in  the  forcing  function.  In  large  problems  it  is  very  important 
from  a  cost  standpoint  to  take  advantage  of  any  symmetry  that  exists. 
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IV  COUPLED  EQUATIONS 


In  this  section  it  will  be  shown  how  the  finite  element  equations  of  Section  II  and 
the  acoustic  equations  of  Section  III  can  be  combined  in  order  to  solve  the  coupled 
structure-acoustic  problem.  This  approach  h&s  been  implemented  at  NOSC  and  has 
been  successfully  applied  to  a  variety  of  structural  acoustic  problems. 


Making  use  of  the  piecewise  constant  approximation  to  p  in  equation  (9),  one 
obtains 


Fm  =  ]£pn  L  4m  ' 

n=l  JSn 

(26) 

Equation  (26) 

can  be  written  in  the  matrix  form 

Fr  =  -CDP, 

(27) 

where 

Cmn  =  i-  1  4>mrtdS 

(28) 

and 

D  =  diag(5j, . .  -,Sn)- 

(29) 

It  is  convenient  to  choose  the  finite  element  break-up  so  that  the  portions  of  the 
elements  lying  on  the  wet  surface  5  coincide  with  the  CHIEF  regions  Sj,...,5/v. 
If  this  is  the  case,  then  the  matrix  C  can  be  computed  within  the  finite  element 
program.  The  n-th  column  of  C  corresponds  to  a  load  vector  in  which  a  constant 
pressure  — 1/5„  is  applied  on  Sn  with  zero  pressure  being  applied  to  the  remainder  of 
the  surface. 

Since  CHIEF  and  the  finite  element  program  use  different  orders  of  interpolation, 
it  is  not  possible  to  fully  impose  the  constraint  for  continuity  of  normal  velocity.  One 
approximate  method  for  imposing  this  constraint  is  to  equate  the  normal  velocity  vn 
in  the  acoustic  equations  to  the  average  of  the  normal  component  of  the  finite  element 
velocity  over  Sn,  i.e., 

Vn  =  4"  /  Mu(x)  •  n{x)\ds.  (30) 

Ofi  JSn 
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Substitution  of  the  interpolation  relation  (1)  into  equation  (30)  gives 


v„  =  [  d>m{x)-n(x)ds) 

Jsn  J 

C^m.  (31) 

m=  1 

Equation  (31)  can  be  written  in  the  matrix  form 

V  =  iuCTU.  (32) 

Solution  of  equation  (7)  for  U  gives 

U  =  (—u2M  +  K)~XFR  +  (~u2M  +  K)-'Fd.  (33) 

Substitution  of  (27)  into  (33)  yields 

U  =  —{—u2M  +  K)~'CDP  +  (—u>2M  +  K)-'Fd.  (34) 

Combination  of  equations  (32)  and  (34)  results  in 

V  =  -iu>CT(-u>2M  +  KY'CDP  +  iu>CT  (—u>2  M  +  K)~XFD.  (35) 
Substitution  of  (35)  into  (15)  gives 


AP  =  -iwflCr(-u5M  +  K)~XCDP  +  iuBCT{-u>2M  +  K)~x  F°  +  Ptne,  (36) 

or  equivalently, 

[A  +  iuBCT{-u2M  +  K)~XCD)P  =  iuBCT(-u2M  +  K  )~XFD  +  Pinc.  (37) 

If  FD  and  P,ne  are  prescribed,  then  P  can  be  obtained  by  solving  the  (possibly  over- 
determined)  system  of  equations  in  (37).  Once  P  is  obtained,  V  can  be  obtained  from 
(35). 

In  the  CHIEF  program  the  system  of  equations  (37)  is  solved  in  the  least  squares 
sense  using  Householder  transformations10.  To  solve  an  overdetermined  system  of 
equations 

MX  =  b.  (38) 

10P.  Businger  and  G.H.  Golub,  Linear  Least  Squares  Solutions  by  Householder  Transforma- 
lions, Numer.  Math.  7,  269-276  (1965). 
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by  this  method,  a  unitary  matrix  Q  is  constructed  so  that 

QM  =  f  U0  )  ,  (39) 

where  \J  is  a  square  upper  triangular  matrix.  The  matrix  Q  is  a  product  of  House¬ 
holder  transformations.  The  least  square  solution  to  (38)  is  given  by 

X  =  (U~\0)Qb.  (40) 


If  equation  (37)  is  to  be  solved  at  a  number  of  frequencies,  then  it  is  more  eco¬ 
nomical  to  obtain  the  inverse  (—  u>3Af  -}-  K)~l  in  terms  of  normal  modes.  The  normal 
modes  Em  (m  =  1, . . . ,  L)  are  solutions  of  the  generalized  eigenproblem 

(~<M  +  E)Em=0,  (41) 

where  w u;£  are  the  associated  eigenvalues.  It  can  be  shown  that  the  normal 
modes  are  M-orthogonal,  i.e., 

E^MEn-0,  m^n.  (42) 

The  modes  are  usually  normalized  so  that 

ETmMEm  =  1,  m  =  1 . L.  (43) 

Let  E  be  the  L  x  L  matrix  whose  columns  are  the  normal  modes  Ei, . . . ,  El.  Then 
the  relations  (42)  and  (43)  can  be  expressed  in  the  single  matrix  relation 


ETME 

=  I. 

(44) 

It  follows  from  (41)  and  (44)  that 

etke 

=  n, 

(45) 

where  Cl  is  the  diagonal  matrix 

( 

0  ) 

n  = 

■ 

(46) 

i 0 

“1 ) 
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If  the  matrix  T  is  defined  by 


then  it  follows  that 
If  W  is  defined  by 
then  W  satisfies 


T  =  (-w2A/  + AT1, 


(— u;2A/  +  K)T  =  /. 
T  =  EW, 


(47) 

(48) 

(49) 


(-u>2A/  +  K)EW  =  /. 

Multiplying  (50)  by  ET  and  making  use  of  equations  (44)  and  (45),  yields 

(-lo7I  +  SI)W  =  Et. 

It  follows  from  (51)  that 

W  =  (~U2I +  Sl)-'ET, 

and  hence  that 

(-u >2m  +  K)~l  =  t  =  £(-u;2/  +  n)-'£;r. 

Since  —  u>2/  +  fl  is  diagonal,  its  inverse  at  each  frequency  is  simple  to  compute, 
stitution  of  (53)  into  (37)  gives 


(50) 

(51) 

(52) 

(53) 
Sub- 


[A  +  iu>BCTE{-u*I  +  ?l)-'ETCD)P  =  iuBC7 E{—u2 1  +  nyxETFD  +  Pmc.  (54) 
If  X  is  defined  by 

X  =  ErC ,  (55) 

then  equation  (56)  can  be  written  as 

[A  +  iuBXt(-u2I  +  ny'XD}P  =  iuBXT(~uPl  +  ny'ETFD  +  Pinc.  (56) 
Given  F°  and  Pinc ,  equation  (54)  can  be  solved  for  P. 


The  solution  of  the  coupled  problem  can  be  summarized  as  follows: 


1.  Compute  the  matrices  M ,  /f,  and  C  using  the  finite  element  program. 
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2.  Calculate  the  matrix  E  of  normal  modes  and  the  diagonal  matrix  17  of  charac¬ 
teristic  frequencies  by  solving  the  generalized  eigenproblem 


(-ulM  +  K)Em=  0,  m  =  l, . . .  ,L. 

3.  Generate  the  matrix  X  =  ETC. 

4.  Compute  the  surface  matrices  A,  B  and  the  far-field  matrices  AFF,  BFF  using 
CHIEF. 

5.  Form  and  solve  the  system  of  equations  (56)  for  P  at  each  frequency. 

6.  Compute  the  normal  velocity  V  at  each  frequency  using  equation  (35). 

7.  Solve  for  the  scattered  far-field  pressure  vector  PFF  at  each  frequency  using 
equation  (23),  i.e., 

PFF  =  AffP  +  BffV. 
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V  SCATTERING  PROBLEMS 


In  a  pure  acoustical  scattering  problem,  the  mechanical  or  electromechanical  load 
vector  FD  is  zero  so  that  the  forced  vibrations  of  the  structure  are  caused  entirely  by 
an  incident  pressure  wave  p,„c.  In  most  cases  p,ne  is  taken  to  be  a  plane  wave  arriving 
from  a  particular  direction,  but  the  formulation  used  in  this  report  can  be  used  to 
calculate  the  scattering  due  to  any  incident  field 

The  fundamental  quantity  of  interest  is  the  scattered  pressure  p„  which  is  defined 
as  the  difference  between  the  total  pressure  p  and  the  incident  pressure  p,nc.  It  is 
usually  desired  to  know  the  scattered  pressure  on  the  surface  and  in  the  far-field  of 
the  object  ao  well  as  the  velocity  or  displacement  on  the  surface  of  the  object. 

The  basic  computation  procedure  follows  the  eight  steps  outlined  in  Section  IV 
with  FD  set  equal  10  zero.  This  procedure  produces  the  surface  pressures,  forces, 
velocities,  and  displacements  as  well  as  the  scattered  pressures  exterior  to  the  object. 
This  method  accounts  for  the  complete  interaction  of  the  incident  acoustic  field  with 
the  elastic  structure  to  the  extent  that  the  structure  is  accurately  represented  by  the 
finite  element  model  and  that  the  acoustic  field  is  accurately  represented  by  CHIEF. 

The  method  can  be  easily  adapted  to  the  simpler  problem  of  scattering  from  a 
rigid  object  ( V  =  0).  In  this  case,  it  follows  from  equation  (15)  that  P  =  A~xPinc. 
Thus,  the  scattered  pressure  field  can  be  computed  as  in  Section  IV  with  P  =  A~lP,nc 
and  V  =  0. 

Typically,  one  normalizes  the  scattered  pressure  by  the  incident  pressure  and 
calculates  either  the  target  strength,  the  form  function,  or  the  scattering  cross  section. 
The  target  strength  is  defined  by 

Target  Strength  ( TS )  *  10  log,,  .  (57) 

where  the  scattered  intensity  I,  is  the  far-field  intensity  of  the  spherically  spreading 
scattered  wave  normalized  back  to  one  meter  (or  one  yard),  i.e., 
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I,  =  lim  R2I{R) 

R—  oo 


=  lim 


R2\P.(RW 


H—oc  pc 

pc 

and,  for  plane-wave  incidence,  the  incident  intensity  7,nc  is  given  by 

\Ptnc{0,nc,<t>tnc)\2 


7,'nc  — 


Combining  (58)  and  (59)  with  (57),  gives 


pc 


(58) 


(59) 


TS  =  20  log, o 


l/’f'tM)! 


\Pinc(0inci4>inc)\, 


(60) 


The  quantity  TS  in  equations  (57)  and  (60)  is  often  referred  to  as  the  bistatic  target 
strength  whereas  the  particular  case  where  ( 0 ,  <f)  =  (#,n£,  <j>inc)  is  called  the  monostatic 
target  stiength. 

Sometimes  a  linear  rather  than  a  logarithmic  measure  of  scattering  strength  is 
desired.  One  such  measure  is  the  scattering  cross  section  a  which  is  related  to  target 
strength  by 

0,„c,dw) 


or  equivalently 


TS  =  101og10 


a{0,<t>,0,nc,4>,nc)  =  4n 


4  n 


(61) 


(62) 


|  P inci^tncy  fiinc)  | 

Another  commonly  used  linear  scattering  measure  is  the  form  function  f  defined  by 


|/(0,<Mm<M<&nc)|  =  J  — 


p.ff(o^) 


Pinc{@inci  fitnc) 


(63) 
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where  ag  is  the  geometric  (high-frequency)  scattering  cross  section  of  the  rigid  object. 
Thus,  for  a  sphere 


\f{9,<t>,9ine, 


«ncj  ^inc) 


(64) 


All  of  the  above  expressions  are  valid  for  the  general  three-dimensional  bistatic  case 
where  (0,nc,  &nc)  is  an  arbitrary  incidence  direction  and  (9,<f>)  is  an  arbitrary  direction 
of  observation.  Frequently  it  is  of  interest  to  compute  the  so-called  monostatic  (back- 
scattered)  target  strength  where  6  =  6inc  and  </>  =  4>{nc. 

As  an  example  problem,  consider  the  scattering  from  a  hollow  steel  shell  immersed 
in  water  and  having  the  following  parameter  values: 


•  Mean  Radius,  a  =  5  m 

•  Thickness,  h  =  0.15  m 

•  Young’s  Modulus,  E  =  2.07  x  10n  N/m2 

•  Poisson’s  Ratio,  v  —  0.3 

•  Shell  Density,  p,  =  7669  Kg/m 3 

•  Fluid  Density,  pw  =  1000  Kg/m3 

•  Fluid  Sound  Speed,  c*,  =  1524  m/ sec 

•  Structural  Loss  Factor,  tj  =  0,  0.01,  0.05,  0.1. 


An  analytical  solution  to  this  problem  is  given  in  Junger  and  Feit11.  This  solution 
can  be  expressed  in  the  present  notation  as  follows 


l/(M)l  =  s 


E 

n=0 


(2 n  -f  l)Pn(cos  9) 

hn]'(ka) 
Jn(^)  ~ 


Pw^u 


(kayh^\ka)(Zn  +  z„) 


(65) 


nJunger  and  Feit,  Sound,  Structures,  and  Their  Interaction,  Second  Edition,  MIT  Press, 
1986,  p.  354. 
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where  the  structural  modal  impedance  Zn  is  given  by 


_  .  k  [n»  -  (nj,‘))»|[n»  -  (ngy; 

"  ^onJ-(nJ  +  n-l+i/)(l+^)' 
and  the  acoustical  modal  impedance  zn  is  given  by 

h^(ka) 

z.  -  v-c. 

Here  (Q^)2  and  (flj,2))2  are  the  roots  of  the  quadratic  equation  in  ft2 

H4  -  [1  +  3i/  +  An  -  02(l  -  i/  -  A2  -  j/An]ft2 

+(An  -  2)(1  -  *2)  +  P\vl  -  4A2  +  An(5  -  t/2)  -  2(1  -  v>)) 

=  0, 

and 

An  =  n(n  +  l) 

P1  =  /i2/(12a2) 

=  JE/M 1  -!/>)] 

fl  =  uja/cp. 


(66) 


(67) 


(68) 


A  description  of  the  numerical  model  set-up  for  this  problem  is  given  below. 

1.  A  spherical  coordinate  system  was  used  having  its  center  at  the  center  of  the 
sphere  and  its  polar  axis  in  the  direction  of  the  incident  wave.  Thus,  all  fields 
are  axisymmetric  in  the  ^-direction. 

2.  Forty  quadratic  axisymmetric  thick  shell  elements  were  used  in  the  finite  element 
(MARTS AM)  calculations. 

3.  The  CHIEF  computations  were  made  with  10  rotational  symmetry  blocks  hav¬ 
ing  one  element  per  block  in  the  (^-direction  and  40  elements  per  block  in  the 
^-direction.  The  Gaussian  quadrature  order  was  16  in  the  (^direction  and  four 
in  the  0-direction. 
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4.  Since  the  form  function  varies  rapidly  with  frequency  when  the  structural  damp¬ 
ing  is  low,  a  very  fine  resolution  in  frequency  was  used.  The  increment  in  ka 
was  0.01  as  ka  varied  from  1.5  to  3.09. 


Figure  2  shows  a  comparison  of  the  numerical  and  the  series  solution  for  the  form 
function  in  the  back-scatter  direction  versus  ka  when  there  is  no  structural  damping 
(r)  =  0).  The  accuracy  is  quite  good.  Figure  3  shows  a  similar  comparison  in  the 
forward  scatter  direction.  Figures  4  through  6  show  how  the  frequency  response  of 
the  back-scatter  form  function  smooths  out  as  the  damping  coefficient  r)  takes  on  the 
values  0.01,  0.05,  0.1.  Damping  is  included  in  the  models  by  making  the  Young’s 
Modulus  of  the  shell  complex,  i.e.,  E  =  Er(  1  -f  »?).  Since  the  shell  in  this  problem 
is  composed  of  a  single  material,  the  finite  element  stiffness  matrix  is  proportional 
to  E.  Therefore,  the  eigenvectors  are  independent  of  r)  and  the  eigenvalues  are  given 
by  -f  irj)  where  is  an  eigenvalue  for  the  case  rj  =  0.  Figure  7  is  like  figure 

6  except  for  the  forward  scatter  direction.  Finally,  figure  8  shows  how  the  target 
strength  varies  with  angle  for  ka  =  2.37  and  rj  =  0. 
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Figure  6.  Form  function  in  the  back-scatter  direction  (jj  =  0.1). 
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Figure  7.  Form  function  in  the  forward-scatter  direction  (77  =  0.1).  j 
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Figure  8.  Target  strength  versus  angle  (lea  =  2.37,  tj  =  0). 
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VI  RADIATION  PROBLEMS 


In  a  structural  radiation  problem  the  incident  pressure  pinc  is  taken  to  be  zero  and 
the  structure  is  forced  to  vibrate  by  a  prescribed  mechanical  or  electromechanical  load 
vector  FD .  Again  the  basic  computational  steps  outlined  in  Section  IV  are  followed 
with  P,ne  =  0.  Since  Pine  =  0,  the  formulas  for  scattered  pressure  now  apply  to  the 
total  pressure  p  as  well. 


To  assess  the  directional  response  of  a  radiating  structure  one  usually  plots  the 
quantity 


K{9,4>)  =  10  log  10 


PFF(0 .  4>) 


pFF{0o,4>o) 


where  (0O,  4>o)  is  some  prescribed  reference  direction  which  is  usually  taken  to  be  the 
direction  of  maximum  response.  The  quantity  K(0,  <f>)  is  called  the  directivity  pattern 
of  the  radiator. 


The  source  level  (SL)  of  the  radiator  is  defined  by 


SL 


10  log  i0 


lOlogjo 


lim 

R—oo  1  fi  Pa 

PFF(0o,to)  2 
10-6  Pa  ' 


(69) 


The  total  radiated  power  W  is  given  by 


where  /  is  the  intensity  and  Sr  is  a  spherical  surface  of  radius  R  which  completely 
surrounds  the  radiator.  The  directivity  factor  D  is  defined  by 


I(R,0o,<t>o) 

TTmU75 


(71) 
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In  the  far-field  we  have 


2 


I(R-.  Go,  A>) 


\p(RJo,4>o)\2 

Pw^ui 


pFF(0o,  4>o) 


(72) 


Combination  of  equations  (69)-(71)  gives 

SL  =  10  log, 0  (—  ^  p^)  +  10  log, 0  tV  +  Dl,  (73) 

where  DI  —  101og10Z).  The  quantity  Di  is  called  the  directivity  index.  If  pwcw  is 
taken  to  be  1.5  x  10fl  which  is  a  typical  value  for  water,  then 


_ Pu/Cyj 

4tt  *  10-12  Pa 


170.8, 


and  hence  equation  (73)  becomes 


SL  =  170.8  +  10  log  l0W  +  DI. 


(74) 


The  source  level  can  be  computed  using  equation  (69)  and  the  total  radiated  power 
can  be  computed  from 

=  itSmMPmV-J,  (75) 

m  =  l 

where  pm  and  vm  are  the  rms  pressures  and  normal  velocities  on  the  m-th  CHIEF 
subdivision.  Therefore,  the  directivity  index  DI  can  be  computed  using  equation 
(74),  i.e., 

DI  =  SL  -  10  log10  W  -  1 70.8  .  (76) 


As  an  example  problem  consider  the  same  spherical  shell  that  was  used  in  the 
example  scattering  problem  of  Section  V.  This  time  the  incident  field  Pine  will  be  set 
to  zero  and  the  shell  will  be  driven  by  a  constant  unit  pressure  applied  over  a  polar 
cap  having  a  half-angular  width  of  4.5  degrees.  An  analytic  series  solution  for  the 
far-field  pressure  is  given  by 


pce~ikR  ^  tn[Pn_i{cos0o)  -  Pn+,(cos0o)]  D 
"  2  kR  to  h?'(ka)(Zn  +  zn)  "(  (  } 
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where  htf\  Zn  and  zn  are  the  same  as  in  the  scattering  problem  and  0O  is  the  angular 
extent  of  the  driven  spherical  cap. 

Figure  9  shows  a  plot  of  source  level  versus  ka  (ij  =  0)  for  both  the  analytic  series 
solution  and  the  numerical  solution. 

Again  the  agreement  is  quite  good.  Figure  10  shows  the  same  comparison  for 
the  case  of  high  structural  damping  (ij  =  0.1).  Finally,  figures  11  and  12  show  polar 
plots  of  the  normalized  far-field  pressure  versus  angle  for  a  ka  of  2.37  and  a  ka  of  3.0. 
Notice  that  the  shape  of  the  patterns  changes  dramatically  with  frequency  due  to  the 
resonant  structure  of  the  shell. 
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Figure  11.  Normalized  far-field  pressure  versus  angle  (iba  =  2.37). 
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Figure  12.  Normalized  far-field  pressure  versus  angle  ( kn  —  3.0). 


VII  SUMMARY 


A  general  method  for  solving  acoustic  radiation  and  scattering  problems  in  which 
there  is  significant  interaction  and  effect  of  the  fluid  medium  on  the  motion  of  an 
elastic  object  has  been  developed  and  described  mathematically.  The  method  can  be 
described  in  three  major  parts: 


1.  A  finite  element  model  (MARTSAM),  which  was  developed  to  solve  complicated 
structural  problems  based  on  the  fundamental  equations  of  linear  elasticity; 

2.  An  acoustic  radiation  and  scattering  method  (CHIEF),  which  is  based  on  an 
integral  representation  of  the  solution  of  the  Helmholtz  wave  equation  for  a 
homogeneous  fluid; 

3.  Additional  programmed  logic  to  couple  these  programs  together  in  an  efficient 
computational  scheme  to  calculate  the  principal  surface  and  field  quantities  of 
interest  in  a  radiation  or  scattering  problem. 


The  method  is  capable  of  dealing  with  elastic  objects  which  can  be  of  relatively 
arbitrary  shape  and  can  have  complex  internal  elastic  structure.  Structural  damping 
can  be  included  by  assigning  a  (homogeneous)  damping  factor  to  the  elastic  ma¬ 
terial.  The  forcing  function  can  be  an  arbitrary  incident  wave  or  a  mechanical  or 
electromechanical  force  or  a  combination  of  forces. 

A  major  emphasis  of  this  report  is  to  cast  the  governing  equations  and  relation¬ 
ships  in  a  computationally  efficient  form  to  allow  the  spatial  and  spectral  response 
to  be  calculated  with  high  resolution.  (Nevertheless,  the  fact  that  the  equations  are 
discretized  in  terms  of  numerous  areas  which  are  small  in  comparison  with  the  wave¬ 
length  implies  that  the  computational  time  and  cost  is  a  steep  function  of  frequency, 
and  the  method  is  fundamentally  a  low-frequency  method.)  Other  factors  which  have 
been  carefully  considered  in  developing  the  computational  algorithms  are  the  poten¬ 
tial  symmetries  inherent  in  the  object  and/or  the  forcing  terms,  and  the  ability  to 
save  major  intermediate  results  for  recalculation  of  similar  problems  with  different 
parameters. 
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Surface  pressures,  displacements,  and  velocities  as  well  as  field  pressures  can  be 
determined.  For  scattering  problems,  the  output  can  also  consist  of  bistatic  or  mono¬ 
static  target  strengths,  scattering  cross  sections,  or  form  functions.  For  radiation 
problems,  the  source  level,  directivity  patterns,  directivity  index,  and  the  total  ra¬ 
diated  power  can  all  be  calculated,  as  well  as  intermediate  results  such  as  radiation 
impedance. 

The  accuracy  of  the  method  has  been  illustrated  by  calculating  both  radiation 
and  scattering  from  a  submerged  spherical  shell  over  the  frequency  range  where  the 
size  of  the  object  is  comparable  with  the  acoustic  wavelength  and  where  significant 
interaction  and  resonant  effects  are  exhibited.  The  numerical  results  from  this  method 
compare  favorably  with  an  analytic  (series)  solution  for  this  canonical  problem. 

The  method  was  developed,  of  course,  to  obtain  accurate  and  high  resolution 
results  in  much  more  general  or  complicated  problems  which  are  not  analytically 
tractable.  In  fact,  it  has  been  applied  with  success  to  complex  transducer  array 
problems  and  to  determining  the  target  strength  of  realistic  submarine-like  structures 
at  low  frequencies.  Results  of  these  practical  Navy  problems  of  interest  are  not 
contained  in  this  report  but  will  be  reported  separately  in  appropriate  reports  and 
publications. 


